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Two, replica symmetry breaking specific, quantities of the Ising spin glass — the break- 
point x\ of the order parameter function and the Almeida- Thouless line — are calculated 
in six dimensions (the upper critical dimension of the replicated held theory used) , and also 
below and above it. The results comfirm that replica symmetry breaking does exist below 
d = 6, and also the tendency of its escalation for decreasing dimension continues. As a new 
feature, x\ has a nonzero and universal value for d < 6 at criticality. Near six dimensions 
we have X\ c = 3 (6 — d) + 0[(6 — d) 2 ]. A method to expand a generic theory with replica 
equivalence around the replica symmetric one is also demonstrated. 

PACS numbers: 75.10.Nr, 05.10.Cc 



I. INTRODUCTION 

Frustration in disordered systems gives rise to a complex equilibrium state with a nontrivial 
breaking of ergodicity (see for a review and important reprints of the field). In the mean 
field version of the Ising spin glass 0], the decomposition of the Gibbs state into ultrametrically 
organized pure states is (mathematically) encoded in the replica symmetry broken (RSB) solution 
of the replicated system [1]. This solution has characteristics — such as the order parameter 
function q(x), and the spin glass transition in nonzero external magnetic field along the so called 
Almeida-Thouless (AT) line — which fully distinguish it from the much simpler replica symmetric 
(RS) case. This RS solution is unstable in the mean field glassy phase jj. 

From the physical point of view, RSB implies the presence of violations of nontrivial fluctuation- 
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dissipation relations at off-equilibrium (during aging), while the off-equilibrium fluctuation- 
dissipation relations would be trivial in the RS case: in particular no aging of the response function 
is expected then, in variance with the experimental evidence in three dimensions at zero magnetic 
field. It is a very important task to determine the dimensional regime where the low temper- 
ature phase with aging response function survives. Evidently, there is no glassy phase in the 
one-dimensional system, whereas there is an ample numerical evidence against any transition in 
the two-dimensional case too. Generally speaking, we expect that the transitions disappear at 
the corresponding lower critical dimensions, i.e at d^ G in zero magnetic field, and at dg G in the 
presence of a magnetic field. We cannot say a priori if these two lower critical dimensions are the 
same: in the case of an Ising ferromagnet with a random magnetic field, for instance, it is well 
known that dP IF = 1, whereas dj F = 0. The situation in spin glasses is quite unclear: the different 
structure of the low momentum singularities in zero and nonzero magnetic field 4] suggest that 
dfg G < dg G , while the arguments based on domain wall energies give d° SG = d h SG = 2.5 [5]. The 
existence of a low temperature phase with aging response function should be ultimately decided 
by investigating the structure of infrared divergences in the perturbative expansion, and by the 
analysis of nonperturbative contributions. This task goes by far beyond the goals of the present 
paper. We aim to study in details the properties of the low temperature phase near the critical 
temperature, and around the upper critical dimension (i.e. six) where the critical exponents at 
zero magnetic field become nontrivial. Our study also aims to correct some recent claims on the 
nonexistence of a RSB phase below six dimensions that are due to an incorrect analysis of the 
consequences of some renormalization group equations [f| . 

The mean field Ising spin glass, at least when studied with the replica trick, can be consid- 
ered as the infinite-dimensional limit of the replica field theory representing the d-dimensional 
short ranged model defined on a hypercubic lattice [7]. The study of this replica field theory for 
decreasing dimensionalities seems to be a good strategy for reaching a full understanding of the 
three-dimensional Ising spin glass. 

n 

This project has had by now a long history whose first period was summarized in Ref. [8J. It 
turns out from these studies that the RS glassy phase is notoriously unstable even down to d < 6, 
with a persistently escalating RSB phase (see, for instance, Fig. 1 of Ref. [9(). A scaling picture 
was proposed in [10j for helping to understand one-loop calculations in the (zero external magnetic 
field) RSB phase. Some of the results of this reference are reproduced and/or revised in the present 
paper, especially the behaviour of the breakpoint x\ of q{x) around six dimensions. The AT line 
was first found in Ref. [ll] for the range 6 < d < 8, whereas it was followed up from mean field 



(d = oo) to d < 6 (and also for nonzero replica number n) in 

Nevertheless, the RS spin glass phase has remained an alternative due to the so called droplet 
model 
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141 ] . This theory predicts a unique Gibbs state (apart from spin inversion) for T < T c — 
that is why the replicated theory is RS — which is massless, and the glassy phase is unstable for 
any infinitesimal magnetic field, i.e. there is no AT line. A schematic picture of the two scenarios 
on the temperature-magnetic field plane is presented in Fig. [TJ The phase boundary lies along the 
temperature axis in the droplet case, a zero-temperature fixed point governing its behaviour; the 
analogous attractive — and also zero-temperature — fixed point for the RSB scenario is shifted to 
a nonzero external field h c . The other end of the phase boundary is, in both cases, the zero-field 
critical fixed point at T c . Since the symmetry of the transition line — namely, an RS state with 
nonzero order parameter q, which is massless in the so called replicon sector, while massive in the 
longitudinal one — is the same (notwithstanding the fact that the AT line proceeds in nonzero 
magnetic field), the two renormalization group (RG) pictures can be studied in a common field 

fis| . The vicinity 
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theory. This is the generic replica symmetric field theory elaborated in Refs. 
of the (hypothetical) zero temperature fixed point can be studied in this field theory by assuming 
a hard (practically infinite) longitudinal mass, th us p rojecting the theory into the replicon sector. 
This was done decades ago by Bray and Roberts [16(, who found a stable Gaussian fixed point for 
d > 6, whereas it was impossible to find any physically relevant and stable fixed point for d < 6. 
This was later interpreted [jjj] as a sign that the AT line disappears below six dimensions, and 
the droplet scenario takes over. This is, however, a faulty argument, since — as we have explained 
above — the RG equations (those for instance of Ref. [16() are not specific to the low temperature 
behaviour of the AT line. An effort to understand the crossover from the zero-field critical fixed 
point to the zero-temperature one was made in Ref. [l7], where the whole set of RG equations was 
derived in a first order perturbative renormalization. (The Bray-Roberts equations are naturally 
included there.) The runaway flows found were discussed in details in [7], and it was argued in 
this reference that the RG scheme used could not be expected to detect a zero-temperature fixed 
point in epsilon expansion. But again, the lack of a fixed point with infinite longitudinal mass in 
the RG equations valid around the critical point is not specific to spin glasses, and this property 
cannot distinguish between the two rival spin glass theories. 

In a recent paper Moore and Bray suggest a proof that RSB disappears when six dimensions 
is approached form above. They take the d — > 6 + limit of known first order results, using RG 
arguments, for x\ (the breakpoint of the order parameter function) and the AT line, and find both 
going to zero. We reproduce their results in a more complete RG scheme in Sec. HH and show what 
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FIG. 1: Schematic phase diagrams for a d-dimensional Ising spin glass in the temperature-magnetic field 
plane. There is an RSB glassy phase in (a) bordered by the AT line. On the other hand, the glassy phase is 
RS in (b), and lies in the zero-field subspace. Both the AT line and the zero-field glassy phase are represented 
by the same generic replica symmetric field theory with massive longitudinal and massless replicon modes. 

h 2 



RS 




is the fundamental flaw in their argument. At this point, the reader is advised to jump to Fig. EKb) 
in Sec. I VII where X\ is plotted against dimension along with the so called scaling variable, which is 
effectively the relative error of the approximation. The breakpoint x\ is monotonically increasing 
for decreasing dimension as long as the scaling variable is small. This is the range where the 
approximation is valid! However at around d ~ 6.1, the scaling variable starts to steeply increase 
(and actually goes to 1 for d — > 6), simultanously x\ suddenly changes its behaviour, and falls to 



but it 



zero: this is the effect (and a similar scenario for the AT line) that has been found in 
must be clear that these results fall outside the range of validity of the approximate RG equations. 
As a matter of fact, x\ can be calculated directly in d = 6 (Sec. IIIip . its value is shown as the 
horizontal line in Fig. EJb): it is visibly an extrapolation of the curve from the range where the 
approximation is good. (In fact, it is an old wisdom of the RG theories that the upper critical 
dimension requires special care.) There is only one case where the arguments of Ref. [6] are correct 
[and interestingly enough, this is admitted there below Eq. (18) of that reference], namely just at 
criticality. But that yields only the trivial results for the d = 6 system: x% is zero for T = T c , and 
the AT line starts at the origin, i.e. at T = T c and h 2 = 0, and does not say anything about the 
disappearence of RSB. 1 



1 Somewhat surprisingly, [f| neglects discussing and even citing Ref. [3| , where the AT line is followed up from mean 
field to d < 6. Subsection IV CI reconsiders and comfirms the existence of an AT line below six dimensions. 
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The outline of the paper is as follows: Section[TI]is devoted to the study of the dimensional regime 
6 < d < 8, although the perturbative results of subsection III Al are extensively used in later sections 
too. In Sec. 1111} the renormalization group ideas are specifically applied to the d = 6 case, simply 
following the lines explained in classical RG textbooks (see, for instance, [l8]). The breakpoint 
x% and the AT line are calculated at the upper critical dimension, both displaying logarithmic 
temperature corrections. A method for expanding a general (except that replica equivalence is 
assumed) RSB theory around the RS one is presented in Sec. IIV1 and applied to the ultrametric 
case. By this method, quantities of the RSB theory, like xi, can be expressed in terms of vertices 
of the RS theory. In the next section, Sec. El we return to our original program, and study the 
case d < 6: generic RG arguments are presented, and the calculation of x\ and the AT line in e- 
expansion is performed. A new feature emerges below six dimensions, namely x\ becomes nonzero 
and universal at criticality. In the last section, Sec. IVH special examples, both for x\ and the AT 
line, are used to conclude that RSB escalates both in the regime above and below six dimensions. 



II. FORMULATION OF THE SPIN GLASS PROBLEM FOR 6 < d < 8 

The simplest replicated field theory corresponding to the Ising spin glass in zero external mag- 
netic field and below d = 8 has two bare parameters defining the model: r (measuring the distance 
from criticality and w (the only bare cubic coupling compatible with the symmetrical — param- 
agnetic — state). Its Lagrangian is 

c = \Y.(\r*+AY. <K - ^ « E' E <«« d) 

P a/3 P1P2P3 a/?7 

where the bare mass fh = fh c — r, and the critical mass has been presented in the literature several 
times in leading order of the loop expansion: 

1 . . 2 1 1 
fh c =-(n-2)w ■ 

p 

In this n(n — l)/2 component field theory the fluctuating fields are symmetric in the replica indices 
with zero diagonals: (j)p^ = c/)p a and (j)" a = 0, a,/3 = l,...,n. [Momentum conservation is 
indicated by the primed summation. The number of the Ising spins becomes infinite in the 
thermodynamic limit, rendering summations to integrals over the continuum of momenta in the 
diagrams of the perturbative expansion. A momentum cutoff A is always understood to block 
ultraviolet divergences, although it can be (and will be) absorbed into the definition of different 
quantities.] The replica number n goes to zero in the spin glass limit. 
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A. Perturbative results 

We are now going to recollect several results for the replica symmetric (RS) spin glass phase 
- see Refs. [j], Q, Q, [3] — which are needed for the following discussion. Due to the severe 
technical difficulties, only one- loop calculations have been accomplished (e = 6 — d and n = 0). 

• RS order parameter q, i.e. the equation of state: 

A 

^ - 1 - 2w 2 r^ 2 - V p2 ~ 2 + -w t- 2 h 2 (2) 

(The last term with the external magnetic field h has been included here for later reference. 
At the moment, it is to be considered as zero.) We can use wq = r in the one-loop diagrams, 
and after rescaling the momentum as p — > pl^fr, two different propagators remain: the 
replicon (p~ 2 ) and the longitudinal [{p 2 + 2) _1 ] ones. To make the formulae for the one-loop 
vertices more transparent, it is useful to introduce a common notation I... for the occuring 
integrals, as is illustrated below: 

A 

_ i^ i rv? d d P i rjt dpp- l+d 

RRLL = N p 4( p 2 + 2)2 ~ J ( 2?r )d p 4( p 2 + 2 )2 ~ d J p 4( p 2 + 2 )2 " 

• The replicon mass: 

T R = 2 mi = -2r + 2wq + Aw 2 t 1+ ^ 2 (4J rll - 3I RRL ) . (3) 

• The basic cubic vertex of the Tr 3 operator: 

wi=w + 2w 3 t^ 2 (-8I RRL + 7I RRR - UIrbll - 8I R lll) . (4) 

• The quartic vertex of 4> a ^ A : 

u 2 = 2Aw a t- 1+ ^ 2 I rrll . (5) 



In fact, this last result is new. Details of the somewhat lengthy calculation of the replicon- type 
quartic vertices will be published later. 
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B. Simple two-parameter renormalization group 

An extensive renormalization group (RG) study of the generic RS glassy phase was published in 
Ref. [17]. When close to the Gaussian fixed point 2 , i.e. w <C 1 and t<1, and only infinitesimally 
breaking the high-temperature (paramagnetic) symmetry of the system, we have the following 
simple two-parameter RG flow-equations: 

'o i i 2 n 4 

w z = —\e\w — 2w , 

10 A < 6 > 

¥ J T - 

Physical quantities take simple scaling forms when, instead of w and r, they are expressed in terms 
of the nonlinear scaling fields w and r defined by: 

~ 2 ii~2 

w = —\e\w , 

(7) 

r = 2r . 

A straightforward calculation provides: 

„r.2\ -i 



w 2 = w 2 I 1 - 2 



.w 

-} 

/ W Z \ ' 

T=r v~ 2 w ■ 

We are now going to compute the quantities q, Tr, w\ and U2 by the RG in terms of w 2 and r. In 
this way, we can get more general results when approaching dimension six from above as compared 
with the perturbative computation: now we may have |e| <C w 2 <^ 1, although the scaling variable 

^2 r |e|/2 mus ^ sma n : 

yj 2 r \ e \/ 2 <g | e | ; even when |e| <C w 2 . 
• The renormalization flow equation for q is 

with t\l = t]r = — ^w 2 in this approximation. It can be solved by using Eqs. ([7]) and ([8]): 

/ ~2\- 1 /6 

q = r 1+ ^q(w 2 r^ 2 ) (l-2^-r) , (10) 



From now on, we redefine the parameters by suitably absorbing the geometrical factor Kd and A: r/A 2 — > t, 
w 2 K d A^ -> «, 2 and ^KJ 1/2 A- 4 -l £ l/ 2 -> h 2 . 



and a comparison with ([2]) makes it possible — after some manipulations — to get the 
leading terms of the scaling function: 

q(x) = -^(l + Cx + ...), with the constant C = 2 1+J ^ r(l + — )r(l - — ) ( + 1 J . 
y x 2 2 \\e\ J 

(11) 

The replicon mass evolves under renormalization as 

£r = (2 - m) ?r = (2 + \ W A t r , (12) 



with the solution 

T R = rf R (w 2 r^ 2 ) (l-2^r J 
Substituting q in Eq. © by r from provides: 



(13) 



T R = -16W 2 T 1+] 4 I RRLL + wr^h 2 . (14) 

Keeping in mind that (|14|) is valid for w 2 ~ w 2 <C |e| and h 2 is zero at the moment, it is 
straightforward to derive the scaling function in Eq. f)13|) : 

t R (x) = c'x + ..., with c" = -2 2 +^r(i + ^)r(i-M). (15) 



|e| 3 \ / |e, 2 



As for w%, we have 

wx=\-^-^tir\wi=\-^+w*\w 1 (16) 
and 

wi = r"x ^ (zZ> 2 ^l-l/ 2 ) (\ - V2 . (17) 

Comparing (JH) and dTTJ) yields 
w 1 (x) = VZ(l + C"x + ...), with C" = 2^r(l + ^)r(l-^) ^-9-|e|V (18) 

Finally, from the flow 

u 2 = ( - 2 - |e| - 27/*) u 2 = (-2 - |e| + ^ 2 ^) u 2 (19) 
follows the scaling form of the most important quartic vertex: 

/ ~2\ 2 /3 

U2 = r-i-J? & 2 (^ 2 r H/ 2 ) ( 1 - 2^ ) . (20) 



From ([5]) and d20j) results [see also (US])] 



u 2 (x) = -^C"x 2 + ... . (21) 
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C. The calculation of x\ and the Almeida Thouless line 

The leading contribution to the breakpoint of the order parameter function q(x) is derived in 
Sec. OVp . and has the simple form [see (|42h and the more general considerations in that section 
about getting x\ on the basis of the generic RS field theory]: 

x\ = — q. 

Inserting (JTUJ) , (fT7|) and (|2"U|) . the scaling equation of x\ follows: 

Xl = xi(w 2 r^^ 2 ), with . . ) = ^ 2 " ' g(. . . ) . 



By the help of Eqs. dTTj) , (HID, PI]) and (USD, we can conclude 

i 1 = 6x2T r(l + ^)r(l - M) w 2 r^ 2 + ... . (22) 
Inverting jSJ), x\ can be expressed by the original bare coupling w: 



Xl ~-^ r W 2 . (23) 

l + 2 v\ 



This equation agrees with Eq. (21) of Ref. |f|. 3 The range of applicability of the above equation: 

w 2 ,r < 1, < |e| < 2 and (most importantly) w 2 r^^ 2 <C |e| . (24) 

If we fix the system's bare coupling w and approach six dimensions, then w 2 — > |e|/2 and x\ ~ 
|e| rl e l/ 2 . This behaviour was interpreted by the authors of Ref. [6] as the sign of the end of RSB 
at six dimensions: a vanishing x\ is consistent with RS. But, as Eq. (|24[) clearly shows, in this 
limit r must go to zero, 4 i.e. the breakpoint disappears at the critical surface in six dimensions — 
a property valid also for d > 6 (but, as we will see later, not for d < 6). 5 In the next section we 
will show that below the critical surface x\ > and has a logarithmic temperature dependence at 
exactly six dimensions. 

We now turn to the problem of the Almeida-Thouless line. The introduction of a magnetic field 
h 2 involves a new nonlinear scaling field h? with 



3 |r(0)| — r in that paper is what we call r here, whereas io(0) = to agrees with our notation for the bare cubic 
coupling. 

4 I — h 

That point has been noticed in Ref. [g], but was completely misinterpreted. We will return to this problem in Sec. 
IVII see the first row of Eq. (|59[1 showing the impossibility of the limit |e| —5- in this approximation. 

5 The multiplicative factor |e| in x\ has its origin in the termination of the definition of the nonlinear scaling field w 
in d = 6. This is a feature of the RS renormalization group, and is not related to the problem of replica symmetry 
breaking. 
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Eq. (I13p remains valid, but the scaling function Tr has now two arguments: x = w 2 r^ e ^ 2 and 
y = h? r~ 2 ~\ e \l 4: . Realizing that the replicon mass starts at one-loop order, the bare parameters in 
(|14p can be replaced by their corresponding nonlinear scaling fields, making it possible to read off 
the scaling function: 

Fii{x,y) = C'x + s/xy; 
see also f|15j) . The vanishing replicon mass defines the AT line, i.e. y = —C \fx providing 

h?= -C'wr 2+ ^ . (25) 
The connection between h 2 and h 2 may be found from the flow equation 

A2 =( 4 + ^T) ft2 = ( 4 + ^ + H ft2 ' (26) 
with the solution [see also ([7]) and (|SJ)]: 

h 2 = h*(l-2?L) . (27) 



|e 

It is useful to display the AT line (j25|) in the original bare parameters by Eqs. (J8]) and ([2 

" 2 = - c '- / 7^Wi r2+ ' 4 - (28 » 

(l + 2f) 

This equation is identical with Eq. (15) of Ref.[6j, and the |e| 4 factor, arising when |e| — > while 
fixing w, led those authors to conclude that the AT line disappears in six dimensions. But, again, 
Eq. (I24D and the discussion below it shows that this limit provides results only on the critical 
surface (r and r zero), and it informs us only about the trivial fact that the AT line starts at the 
origin of the r, h? plain. 

III. AT THE UPPER CRITICAL DIMENSION: d = 6 

As can be seen from the previous section, knowledge about the six dimensional system cannot be 
gained from the RG results in the d <J 6 case. The fundamental reason for that is the impossibility 
to linearize the RG flow equations at exactly an upper critical dimension. Therefore, the scaling 
field w is not defined for d = 6, and we keep w (although r and h? are still meaningful). The RG 
flow ([6]) is now: 

w 2 = —2w , 



10 



2 w z t . 



(29) 
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The connection between r and r becomes [instead of ([Sj)] : 



10 



(30) 



and the scaling variable with zero scaling dimension is now (instead of w 2 r' e '/ 2 ): 

,2 



■ur 



1 — w 2 In r ' 

which can be easily checked by Eq. (|29l) and the nonlinear scaling field property r = 2r. 

A. The calculation of x\ 



The renormalization group flow equations for the three relevant physical quantities q, w\ and 
U2 are as follows: 

1 



q=[2--W 



W\= W Wl , 
U 2 = I -2 + ^-W Z I U2 ■ 



They all have the same form, and their solutions are easily found in scaling form. 



The RS order parameter: 



q = ui3 rq 



w 



1 — w 2 In r / ' 



j(x) = x 



1 + (2 + In 2) x + - x In x + . . . 



(31) 



The scaling function q(x) has been obtained by evaluating ([2]) in d = 6 (in zero magnetic 
field at the moment) and using the connection between r and r in (|30p . 

The cubic vertex ioi in six dimensions: 



W\ = w 1 Wl 



ur 



1 — w 2 In ', 



rt?l = a; 



1 + 



^ + 8 In 2 — 71nn N )x+^xlnx + ... 



(32) 



Eqs. ([!]) and (I30p has been used to get the scaling function. One important remark is 
appropriate here. The term with the logarithm of the replica number, Inn, comes from 
Irrr i n ©j and is a prominent example of the severe infrared divergences caused by the 
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replicon propagator. Similar contributions enter in higher order vertices, such as Irrrr in 
the quartic vertex belonging to the operator Tr<^ 4 . This is a clear indication — beside the 
instability of the replicon mode — that the replica symmetric theory is ill-defined in the 
spin glass limit. In fact, these infrared divergent terms can be resummed when we build up 
the RSB theory on the basis of the RS one, as explained in Sec. IIVI What is gained in this 
resummation, after setting n to zero, is the small mass regime of the RSB solution which 
effectively acts as an infrared cutoff. It must be stressed that without this resummation, the 
theory is infrared divergent in any arbitrarily high dimension. 

As for the quartic vertex 112, its scaling form and the leading term of the scaling function 
are [see P and (l30l)]: 

_i 1 „ / w 2 
U2 = w 3 r U2 



1-t^lnr, (33) 
U2(x) = 6x + . . . . 

By Eqs. (I31D . (132)) and (I33D x\ turns out to be a function of the scaling variable, as it must be: 

U 2 „ / W 2 \ „ . . U2{---) ~, s. 

xi = — q = xi T — with xi{. ..) = — -g ... . 

The leading order of the scaling function is simply xi(x) = Gx + . . . , providing one of our basic 
results 

w 2 \ . 10 



X\ = 6 ( 7T- — + . . . ; ffi,r<l and r = tw 3, d = 6 . (34) 

\l-rmr/ 

It is clear from the above equation that x\ is zero at criticality (r = r = 0), and for fixed w the 
approach to zero is logarithmic: 

x\ = 6 I In rj^ 1 + . . . ; r, r — > and w fixed, d = 6 . 



B. Almeida Thouless line in six dimensions 

The flow equation for the replicon mass is unchanged as compared with the d > 6 case, and is 
given by Eq. f)12|) . The nonlinear scaling field corresponding to the external magnetic field satisfies 
h 2 = 4h 2 , therefore the second variable with zero scaling dimension is h 2 /r 2 . Straightforward 
considerations lead us to 

2 - / w 2 h 2 \ , , 

rR = w -srT R _,_ . 35 

\ 1 — w z In r r z I 
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The evolution of the "bare" magnetic field, i.e. h 2 = (4 + ±w 2 ) h 2 [see ([26])] and ([29]) yield 

h 2 = h 2 w-\. (36) 

Evaluating Eq. (|14p at d = 6, and replacing the bare parameters r and /i 2 by r and /i 2 according 
to (pOj) and ([36]) . respectively, makes it possible to read off the scaling function in leading order: 

T R (x,y) = -(-Ax 4 + y + ...). 
x 

From its zero, the AT line is obtained as follows: 

(2 \ 4 
^75-; — ) +•••; w,r<^.l and r = Tw~^, h 2 =h 2 w^; d = 6 . (37) 
1 — w z In r / 

For a given cubic coupling u;, the magnetic field vs. temperature relationship for the boundary of 
the RS phase when approaching the critical point becomes: 

/i 2 =4r 2 |lnr|- 4 + ... ; r,r — >0 and w fixed, d = 6 . 

IV. FORMULATION OF REPLICA SYMMETRY BREAKING ON THE BASIS OF THE 

GENERIC REPLICA SYMMETRIC THEORY 

The considerations in this section are quite general and, therefore, the paramagnetic system 
(i.e. an RS system with zero order parameter) must be represented — instead of the simple case 
of (prj which is sufficient around d = 6 — by a model which includes all the invariants compatible 



with its higher symmetry 201 ] . The replicated field theory is now defined by the Lagrangian C of 
the symmetrical (high-temperature and zero-field) theory: 

£ -\ E (? 2 + *.) E <K - i^Ti E * E «?« - 2i~jv E 

P ^ ' a/3 P1P2P3 a/37 P1P2P3P4 

a/3-yS a/3 a/37 

- E <««c + **E 

P1P2P3P4P5 a/3-yS/i afi'y 

-3 E <<t&f£f&f& + * E <«CCC) + • • • (38) 

a/375 afi'Yflv 

where the bare mass mi = rh = rh c — r, with r measuring the distance from criticality, has been 
also used in (pQ), and tui = to (momentum conservation is indicated by the primed summations). 
The fifth order invariants with the v bare couplings were also included here. In what follows, we 
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use the same notation for an exact vertex (e.g. 112) and its corresponding bare coupling (^2), the 
bar indicating always a bare quantity. 

As explained in details in Appendix D of Ref. [9] , the generic Legendre-transformed free energy 
can be expanded around the RS spin glass state with order parameter q; see (D.5) of this reference: 



1 TV n 1 

+ 2 



a(3 



+ m 2 ^2(q a -y - q){qp 1 - q) + m 3 (fe/3 - <?)(<? 7 <5 ~ ?) 



wi ^2(q a p - q)(qp-y - q){q ia -q) + w 2 ^(qap - qf + w 3 ^2(q a p - q) 2 (q 

a/37 a fi a fii 



a 7 



+ ... 



1 

24 



ui ^2 (q a p - q)(qfa - q)(q 1 s - q)(qs a -q) + u 2 ^(qa/3 - qf + ■■■ 



aft 



1 

120 



vi ^2 ( q °>P ~ vfvp-y ~ ?)(*y<s ~ <i)(i8n - q){q^ a -q) + 



+ 



(39) 



In zero external field J~(q a p) has the same symmetry as C of Eq. (|38p — which is higher than that 
of a generic RS system — , even when T < T c , and using this symmetry a set of equations can be 



found bolwtvn the exact vertices of thr generic .RS theory (see Rels. j3, Q)- ,n ost clfeetivo 

way to get the required vertex relationships is demanding that invariants incompatible with the 
symmetrical theory, e.g. X^a/3 ?a/3' mus ^ finally disappear from (|39p . In this manner, all the vertices 
of the lower symmetry: m,2, m 3 ; w 2 ,- ■ ■ uig; U5,. . . U23; ■ ■ ■ etc., (see Appendix A of Ref. [9} for the 
classification of cubic and quartic vertices) and, as a bonus, m\ can be expressed in terms of wi, 
iti, u 2 , us, Ui, vi, v 2 , V3, v^, and higher order symmetrical vertices. 6 We than have 

11 1 

mi = -nwi q + -(n 2 ui - 2u 2 ) q 2 + — n{n 2 vi - 2v 2 ) q 3 + • • • , 
2 b 24 

m 2 = -wi q - ~(raii + u s ) q 2 + ^-[5n(3n 2 - 5n + l)vi + 2v 2 - 4nu 3 ] q 3 + . . . , 
3 60 



1 



1 



m 3 



-(ui + 2u A ) q 1 - — [5(5n - 4)ui + 2v 3 + 6nv 4 ] ^ + 



60 



6 A vertex is called symmetrical if it is nonzero in the zero order parameter RS system. 
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and furthermore 

1 2 
w 2 =u 2 q+ ^ nv 2 q + • • • , 

w 3 = u 3 q + ^(3v 2 + nv 3 ) q 2 + . . . , 

w 4 =U4q+ ^(v 3 + 3nv A ) q 2 + . . . , 

w 5 =uiq + ^(5nvi + iv 3 ) q 2 + 

1 2 

^6 = -^v 3 q + 

w 7 = 4g(10ui + 12u 4 ) g 2 + • • • , 

w 8 = 0(g 3 )- 

Of the quartic vertices, only those are listed below which are required up to the order of the present 
calculation: 

3 2 2 

115 = -v 2 q+ u 6 = ~v 3 q + . . . , u 7 = -« 4 q + . . . , 

2 12 
"8 = g«2 9 + • • • j Uio = -t>3 g + • • • , uu = -v 3 q+ 

3 

"14 = 7^4 q + uie = vi q + . . . . 

5 

By exploiting these expressions, the free energy functional in Eq. (|39|) can now be written (omitting 
an additive term depending only on q): 

■JjF{q a fi) = \ M( l ^2 ~ \ W ^2 <?a/3<?/3 7 <? 7 a ~ ^ [(«1 + «lQ + • • • ) X] ^apQp^SQSa 
a/3 aP^S 

2 3 3 / 

+ («2 + g ^2<? + • • • ) 2 9a/3 + («3 + ^V 3 q + • • • ) ^ ^7^7 + 0*4 + g «49 +■■■)(£ ^a/3, 

a/3 a/37 Q /5 

1 r 

(«! + ...) ^ qapqpjq-fSqSnQva + (v 2 + ■ ■ ■) ^ 9a/3<?a 7 <?/3 7 + («3 + • • • ) ^Q^Q-ya^S 
a/3-ySfi a(3f a/3 7 <5 



120 



+ («4 + ...)(E 



a/3 7 a/3 

with the following notations for M and W: 
M = (n — 2)wi 

+ \[{n 2 - 3)ui + u 2 + (n - 1)%] 9 + ^(n 3 - 4) Vl - 2(2n - 8)w 2 + 2(n - l)(n + 4)x5 3 ] </ 2 + • • • , 
3 bU 

W = u>! + ui g + 2q[10«i - 3v 2 - (n - 1)5 3 ] ? 2 + • • • • 
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The tilded vertices M3 = 113 + nu^ and V3 = V3 + nv^ were introduced here; in fact, only these 
combinations will enter the equation of state. 

Stationarity of the free energy functional provides the equation of state: 

1 2 
= Mqq ap - W {q 2 ) aj3 - - [(m + Vl q+...) {q 3 ) a p + (u 2 + -V 2 q + . . . ) q% 

+ {u3 + -v%q+. . . ) {q 2 ) aa q a p] _ go( 5 ( u i+- • • ) (<? 4 W + (^2 + - • • ) [3(<7 2 ) Q/ g ^+^(^7^7+^7907)] 

7 

+ (v3 + ...)[2(q 3 ) aa q a p + 3{q 2 ) aa {q 2 ) a p]} . (40) 

Only replica equivalence was used in the derivation of this equation — (g 2 ) att , for instance, is 
independent of the replica number — , otherwise it is quite general: it provides an RSB solution in 
terms of the RS order parameter q (which measures the distance from criticality now) , and of the 
exact RS vertices. It can equally be used in any regime where some kind of perturbation theory is 
valid. 

We now turn to the case of infinite step, ultrametrically organized RSB. The small parameter 
making possible a perturbative treatment is x±, the breakpoint of the order parameter function: 
it is proportional to q in the SK model and for the field theory above 8 dimensions, to q d / 2 ~ 3 
between 6 and 8 dimensions, whereas it is of order e below 6 dimensions. q(x), the order parameter 
function, has the form: 7 

q(x) = qi [r + x\ 8q(r)], with r = x/x\ and 5q(l) = 0. (41) 

The contributions of the various vertices to Eq. (|40p are listed below. The definition of the 
bilinear expression {...;...} used extensively in that list is as follows: 

{f(r);g(r)} = f(r)g(l) + f(l)g(r)-f(r) f dug(u)-g{r) C duf(u)-rf(r)g(r)- f duf{u)g{u). 

Jr Jr JO 

2{qi -q)r- xiqi (r - ^r 3 ) + 2x\{qx - q) 5q(r) - 2xfqi {r; 5q(r)} + 0{x\ ), 

• u±: 

-(<7i -q) 2 r + xi(gi - q)qi (r - K 3 ) - -x\q\ (-r - ^r 3 + ^r 5 ) + O(xf), 



7 We hope that the ratio r of x to xi introduced here cannot be confused with the temperature-like scaling field of 
previous sections. 
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u 2 : 



\q 2 r - ^q't r A + ^x{q A 5q{r) - x{q{ r l bq{r) + 0(xf), 



,2„3 i 1 _2„2 



U 3 
1 



' / 2 2\ 2 2 1 9/2 



9 9 
2 3 2 3 2 4 

g ) 5g(r) - -x^ r{r; <5g(r)} r= i - -x^ 5q(r) + O^) 



ui: 



^Wl 



^i(?l - qfqi (r - ^r 3 ) + ~xl(q x - q)q\ (~r - ^r 3 + ^r 5 ) 

- — x?g 3 f-r - -r 3 + — r 5 - — r 7 ) + 0(xf ), 
12 m v 2 2 10 14 7 v 17 



• v 2 : 



1/ \ 2 / \ 2 3 H / ^ 3 9 c \ 3 2 / ^ 3\ 

- g wi - q)q r + — Wi - q)qi r - 3^ x ^i + 2 r - ^ r ) + qi ( r - ^ r ) 

- 7^ x iqi r2 ( r ~ l r3 ) + -k(qi ~ qfqr + — Wl - g) 3 r - - g)g? (1 - 2r 2 )<5g(r) 
20 3 10 30 5 

11 3 1 1 



1 

To 



xlql{r;r 2 8q(r)} + 0(x$), 



V3- 



3 2 1 2 

— (gi - g) 2 gi + j^iWi ~ ?)«? - 75^1 + 15 Wi ~ ^ 



1 

H xigi 

20 



2(gi - g)gi - \x iq { - ( qi - g) 2 ] (r - ±r 3 ) + O(xf) 



Inserting the above expressions into Eq. (|40p and demanding that the coefficients of r and r 3 
disappear, xi can be read off with some effort. It is best to give x\ as the zero, /(xi) = 0, of the 
following function: 



( — q)+^{— q )+■■■ 



+ 



l-^g 2 )+... 



x+ 



60 

The leading contribution is the well-known formula 



1 1 / Ul X 

3 6 wi ' 





r 1 1 


x 2 + 


[-9+-J 



x J +. . . . 



U 2 

xi = — q 

Wl 



(42) 
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which is used extensively throughout this paper. As a byproduct, the shift of q\ from the RS order 
parameter is given by 

1 2 

11 - Q = 3 xxq{l + -xx + ...) . (43) 

[To preserve consistency, a sixth order contribution — 2/2 Saff^"/? - should have been included 
in the free energy expansion (|39p . as it enters the constant of f(x) at the third order, i.e. at the 
highest order studied here.] 



V. BELOW SIX DIMENSIONS 
A. The renormalization group: fixed point and nonlinear scaling fields 

In d = 6 — e the Gaussian fixed point becomes unstable, and the zero field spin glass transition 
is governed by the non-trivial one. Here we collect and present the available results for this fixed 
point (in the results for the fixed point below, a generic n is kept, although n = is taken in the 
further parts of the section) : 

wl 2 = w* 2 = e , fh\ = m* = — w* 2 ; see Refs. I21I and [l?! • 

2 — n 4 

Although they will not be used in this paper, the fixed point values of the quartic couplings (which 
- according to our knowledge — have not been published before) are also listed here: 

q Q 

— * *4 — * 1 o *4 — * o a *4 — * *4 

u 1 = -nw , u 2 = 12 w , U3 = — 24 iu , u\ = - w . 

The renormalization flow equations for the bare couplings of the generic RS theory were displayed 
in Ref. Using these equations, a new set of parameters gi — the so called nonlinear scaling 



fields introduced by Wegner 22| — can be defined with the following properties: 8 



gi = at the fixed point for all i. 

An infinitesimally small gi, with all the others kept zero, gives an eigenvector belonging to 
the eigenvalue Aj of the linearized renormalization group equations around the fixed point. 

They satisfy exactly the equations gi = \igi. 



The summary presented in this paragraph about the use of nonlinear scaling fields is quite general, not limited to 
the nontrivial fixed point of the RS replica field theory. 
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The RG flow of an observable y — the order parameter or an irreducible vertex, for instance - 
can be written in terms of the g^s as follows: 

y = fk + ^2hgi + ^2 hj gigj + . . . ) y . (44) 

i ij 

The solution of this equation, i.e. y in terms of the scaling fields is easily found: 

/ v fc/Ai -/ -A 2 Ai -Ai/Ai \ f sr^ fa , faj 

y{9l,92,---)=9i y(52 5i ,---,9i9i , . . . j x exp { ^ j; 9 » + 2^ \. ~1 \ 9i9j + ■ ■ 

ij 1 3 

(45) 



the scaling function y(. . . ) can be determined by perturbative methods. 

In our two-parameter system defined by r and w the two nonzero scaling fields 9 r = g\ and 
9 = 92 (the notations are chosen to keep connection with previous sections) can be found by 
starting with the RG equations © 10 and taking the temperature- like relevant eigenvalue A r and 



the irrelevant one, A^, from Ref. [17[ as 

\ r = - = 2-\e + ... , Xg = -e + .... (46) 

The bare parameters are then straightforwardly expressed by the scaling fields as 

2 *2 . 9 e / 2 
w = w + — 



1-22 1-22 

(47) 

3X-5/3 



t = r [ 1-2 



B. x\ below the upper critical dimension 



For the calculation of x\ to first order in e, the RG study of q, w\ and u<i is required. The 
truncated (one- loop) renormalization group equations Q, (fl~6j) and (p79|) — see also footnote [TUI — 
can be used whenever a) 2 Cl and t<C1. We can solve these truncated equations in a similar way 
as (]45p was derived from the generic equation (|44ji . The scaling exponents and the relations between 
bare and scaling parameters are taken from Eqs. (j46]) and (I47p . respectively. The scaling functions, 
which are always denoted by the "hat" symbol, cannot be determined by the RG equations alone, 
but the perturbative results of Eqs. ((SJ), (HJ) and (0) make it possible to get them to first order in 
e. [The bare values must be replaced by the scaling fields using (|47p. and take into account again 



9 In references 0, H(| an alternative scheme was used with a second relevant scaling field entering after appropriately 
redefining the field theory for getting rid of "tadpole" diagrams. The irreducible vertices are the same in the two 
schemes. 

10 But be careful to replace — |e| with t. 
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footnote [TU1 ] In the following, the results for q, w\ and u 2 are listed in itemized form. The k and 
k 2 = kg quantities defined in (jH]) are also presented for completeness. 



r i+ 2 qtgrz) x 1 - 2 ' 



(48) 



with the scaling function 



w 



q{x)= l + (-ln2 + l)e + 



+ 2 



1 + ( ln2 + ^) e + 



- e )+o 



(49) 



12 X 
A: = 2-| + -r ? 2 = 2-- e + 0(6 2 ) and k~ g = -- + 0(e). (50) 



6 , e. / ^\l/2 

ioi = ra wiygrz ) x (1 — 2 -J , 



(51) 



with the scaling function 



w\{x) 



w 



39 7 

1 + (4 In 2 In n) e + 

v 4 2 ; 



2 1 



56 

i In 2 7 In ra) e + 

3 ; 



e 



(52) 



^=|-^^ = 6 + 0(6 2 ) 



and % = 1 + 0(e) 



-1 / e \ ( 5\ 2 / 3 
u 2 = r 1 u 2 (gr2) x {1-2- ) , 



(53) 



with the scaling function 



u 2 (x) 



w 



*4 



l + 0(e) +12 l + 0(e 



+ o 



(54) 



5 „, 2 . . . 4 



/c = -2 + e- 2n|j = -2+-e + 0(e) and % = - + 0(e) 



The following remarks are appropriate here: Firstly, according to Eqs. (1451) . (I46p and (I50p the 
temperature exponent for the order parameter q is exactly k/X r = (2 — e/2 + ra*/2) f = /3 = 
1 + e/2 + 0(e 2 ); see (08]). The temperature exponents in (@8|), (f5Tj) and (|53|) are correct up to 
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e order. Secondly, the discussion below Eq. (|32p concerning the fully replicon, infrared divergent 
contribution to w% is equally valid for the Inn terms in (I52p , Thirdly, the O(e) corrections in 
the scaling functions q and w\ are unnecessary for the leading order calculation of x\\ they are 
displayed here to demonstrate the general form of the e expansion. The corresponding corrections 
for ii2 are not even available, see (I54D . as they would require a twodoop level calculation. 

The leading contribution in the e expansion for x\ follows from substituting q, w\ and u<2 
from Eqs. (gSJ) -(gSJ, (|5T]) - (|52"]) and (Pj3 ]) - ([5I|) . respectively, into the basic formula in Eq. (@2J). A 
remarkably simple formula reflecting the invariance of x\ under renormalization can be concluded: 

x\ = 6w* 2 x%(gr^) = ?>exi(gr^ , with the scaling function 



Xl ( x ) = l + 0(e) +2 l + 0(e 



+ 



(55) 



C. Almeida-Thouless line for rf < 6 

The external magnetic field h 2 evolves under renormalization according to Eq. (|26p . with |e| 
replaced by — e. The corresponding nonlinear scaling field 53 = h 2 has now the relevant eigenvalue 

h2 2 2 v 

The flow equation for the replicon mass — Eq. (|12p — does not contain explicitly the magnetic 
field, therefore it enters the solution only through the invariant h 2 r~ Sl3 ; see Eqs. (|44p . (|45p and 
(|46p . According to the generic scheme (|45p . we have 

r fi = r( 2 ""*) " f R (g \ h 2 r~^, . . . ) x exp (J A + . . . ) . 

The exponential part can again be calculated in the truncated, onedoop approximation, in the 
usual way, providing (note that Xg = — e + . . . ) 

1 

g\ 3 



1-2 

e 

whereas a comparison with the perturbative result (I14p — after substituting the bare parameters 
by their corresponding nonlinear scaling fields [see Eq. (|47p and also 

^ 2 = ^(l_ 2 f) 1/6 (56) 

which follows from (|26|) ] — gives the scaling function: 

f = u,* 2 {[-4 + 0(e)] + [-24 + 0(e)] (f ) + [l + O^J] (^) + [-2+0( e )] (f ) (^) +• • • } • 

(57) 
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The zero of the scaling function gives the AT-line: 

h? = Aw*r 5p = Aw*r 2+ - , - r e/2 < 1 and 0<e<l. (58) 

e 

For the fixed point system, w = w* implies g = and h 2 = h 2 , r = r. The result in (|58[) is then 
identical with Eq. (18) of Ref. p). 

VI. DISCUSSION OF THE RESULTS AND CONCLUSIONS 

For a thorough analysis of the (^-dependence of x\ while crossing the upper critical dimension, 
we recollect here the one-loop truncated results from previous sections; see Eqs. (|22p . (|34p and (|55p . 
The goodness of these approximations depends on the smallness of the scaling variable, which is 
defined and expressed in terms of the bare parameters w 2 and r as follows: 

2. ^2 r \e\/2 = 2_ w 2 T \e\/2 ^ + 2rfyl-\\*\ d> ^ gee || i 



scaling variable 



w 2 (l - w 2 lnr) 1 = w 2 (l + |w 2 Inw 2 - w 2 lnr) 1 d = 6, see (j30 



5,„*2 



|^ r e/2 = r ^^3 ^_^ f w *2 = f d<6,see(|47j). 

(59) 

This scaling variable is displayed — for a chosen pair of bare values w 2 = 0.005 and r = 0.0001, 
both much smaller than one, as it should be in this approximation — as a function of dimension 
d below [Fig. d2^,)] and above [Fig. ([2b)] 6, where it takes ~ 0.005. x\ is also shown in this figure, 
with the awkward behaviour of approaching zero from both sides of the upper critical dimension 
six, while x\ ~ 0.03 in d = 6. It is clear, however, from the figure that our approximation breaks 
down when approaching d = 6 from either side, as the scaling variable goes to unity in that limit. 
As a matter of fact, it must be stipulated that the scaling variable be at least as good as in d = 6, 
i.e. < 0.005. Therefore, the range of applicability of our approximation (for the given w and r) is 
constrained to d w 5.99 and d > 6.4, respectively. (Note that the chosen w is just the fixed point 
when d = 5.99.) Representative values of x\ in these ranges, together with the six-dimensional 
case, are presented in Tab. HI It can be concluded from this example that x\ keeps on being 
monotonically increasing when lowering dimensions through 6. Nevertheless, a discontinuity of 
xi(w 2 ,t) at d = 6 cannot be excluded. An extrapolation of the data from the range d > 6.4, 
using an exponential and/or a power law fit, provides xi(0. 005, 0.0001) « 0.026 — 0.028, a value 
somewhat lower than the six-dimensional one, 0.030, when considering the scaling variable as a 
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FIG. 2: The scaling variable (left vertical axis) measures the goodness of the approximation, (a): d < 6 and 
(b): d > 6. The dependence of x\ is also shown in both regimes, together with its d = 6 value (horizontal 
lines), w 2 = 0.005 and r = 0.0001 are fixed in this figure. The approximation breaks down when approaching 
d = 6 from both sides. 
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1 — xl for 6<d scaling variable xl,d=tj| 



(a) (b) 

TABLE I: x± around six dimensions shows monotonically increasing behaviour with decreasing dimension- 
ality. The smallness of the scaling variable verifies the approximation. The bare parameters w 2 — 0.005 and 
t = 0.0001 are the same as in Fig. [2] 



d 


X! x 10 2 


scaling variable 


6.8 


0.1287 


0.000308 


6.6 


0.2648 


0.001026 


6.4 


0.5649 


0.003834 


6 


2.9943 


0.004991 


5.99005 


2.9993 


0.004776 


5.99 


3.0000 





5.98995 


3.0006 


-0.004774 



measure of the relative error (it is ~ 0.005 in six dimensions, see Tabled]). A similar extrapolation 
from the d < 6 side, however, does not exist. 

Below six dimensions x\ has only a slight temperature dependence, and it becomes nonzero and 
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FIG. 3: x\ as a function of the reduced temperature r; w 2 — 0.005. (a): d < 6 and (b): d > 6. si is zero at 
criticality, i.e. for r = 0, when d > 6. On the contrary, it is nonzero for d < 6 and has the universal value 
X! =3e + 0(e 2 ) at T c . 



I d=6 d=5.99l 



| — d=6.4 - - d=6.2 d=s1 



0.00005 
tau 



0.00005 
tau 



(a) 



(b) 



universal at criticality: 

Xl = [3 e + o(e 2 )] + Crf+- + ... , d = 6-e; 

C is a nonuniversal, i.e. u>-dependent, amplitude. The typical behaviour for both below and above 
six dimensions is displayed in Fig. [3l the value of the cubic coupling is kept w 2 = 0.005. The 
tendency of an increasing x\ while lowering the dimension is again obvious. The vertical scale 
was magnified in the left subfigure (a) to show the qualitative difference between the 6- and 5.99- 
dimensional curves. 

The critical field along the AT-line, for a given pair of bare parameters w 2 and r, can be analysed 
using results from previous sections. See Eqs. (fT5|) . (|28|) for d > 6, and ([37|) for d = 6. Below six 
dimensions, if we wish to move somewhat away from the fixed point, the zero of the expanded 
equation ([ST|) must be found, providing [instead of 



h? =Aw*r & P (l + 8^r e / 2 ) , 9 - r e / 2 < 1 and 0<e«l. 
e e 

Eqs. (|4"7|) . (|56|) and ([59]) . together with Sf3 = 2 + |e, give us the critical field as 

,*2 _ 6 



tuft =4rf 



U' 2 



V w 2 J V 



it; 

2 
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FIG. 4: Almeida-Thouless line (wh 2 versus r) of the field theoretic model with w 2 — 0.005 for three different 
dimensions. 

"~ — d=B.4 d=6 d=5.99l 



1.5xl0- 10 - 
l.xlO- 10 - 
5.xl0- u - 

-- 



The critical field where RSB sets in as a function of temperature (i.e. the AT line) — or more 
precisely wh 2 as a function of r — is shown in Fig. 2] for three different dimensions at fixed cubic 
coupling w 2 = 0.005. The curve for d = 5.99 (note that the system is at exactly the fixed point 
then) is significantly below the six-dimensional one. It is easy to see that this follows directly from 
the exponent inequality 5(3 — 2 = |e + • • • > 0. To see clearly the behaviour of the critical field 
above and below six dimensions for decreasing d, it is tabulated in Table [U for the system with 
w 2 = 0.005 and r = 0.0001. The last three rows of this table show that the kind of monotonicity 
found above six dimensions is restored below it, i.e. the critical field increases with decreasing 
dimensions. It must be remarked, however, that around the last dimension value d = 5.98995, 
the error 11 of our approximation starts to be the same order of magnitude as the variation of the 
critical field itself. The range where this one-loop perturbative method is applicable below the 
upper critical dimension is certainly very narrow. 

As a conclusion, we can confidently claim that RSB survives below six dimensions in the cubic 
replica field theory representing the Ising spin glass. We focused on two quantities which are 
strongly related to RSB: the breakpoint of the order parameter function x\ and the Almeida- 
Thouless line. A combination of the perturbative one-loop method with a simple two-parameter 
renormalization group (which is correct near the critical fixed point) provided reliable results in 
all the three ranges of dimensionalities, i.e. for d larger, equal, and less than six. The calculations 




1 The relative error can be estimated as being proportional to the square of the scaling variable. 
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TABLE II: Critical field values around six dimensions for w 2 — 0.005 and r = 0.0001. Below the critical 
field replica symmetry is broken. The scaling variable's values are, of course, the same as in Table HI 



d 


wh 2 x I0 10 


scaling variable 


6.8 


0.0827 


0.000308 


6.6 


0.1680 


0.001026 


6.4 


0.3497 


0.003834 


6 


1.9849 


0.004991 


5.99005 


1.7410 


0.004776 


5.99 


1.7419 





5.98995 


1.7421 


-0.004774 



above and below six dimensions are rather different, due to the Gaussian versus nontrivial fixed 
point governing critical behaviour in the two cases. The applied perturbative method makes it 
impossible to approach closely the upper critical dimension: the range of dimensions where the 
approximation is correct for a given system (i.e. for given w and r) is very narrow and close to 6 
when d < 6, whereas it is d > 6.2 when d > 6 (and the farther we are from d = 6, the better the 
approximation). The six-dimensional case needs special care along the way systems at their upper 
critical dimension are commonly studied [la]. The logarithmic temperature dependences obtained 
are quite similar to those in ordinary systems at their upper critical dimension. 

Above six dimensions both x\ and the critical field are monotonically increasing for decreasing 
d, and this tendency persists for d < 6. There seems to be, however, a discontinuity of the critical 
field at d = 6~ : the AT line for d S 6 is significantly below the six-dimensional one, see Fig. 2] and 
Table HH Nevertheless, we can notice that the trend of increasing dominance of RSB for decreasing 
space dimensions persists even below six dimensions. 

As a final remark, we recall that for d < 6, x\ gains the qualitatively new feature of being 
nonzero (and universal!) at criticality. This might suggest a kind of first order transition. That 
this is not the case can be clearly seen by displaying the order parameter function using Eqs. (|4ip 
and (H3D : 

q{x) = qi q(x/x\) , where q\ ~ q ~ t 13 . 

An elaboration of the equation of state along the lines of Sec. IIVI for d < 6 (which is out of the 
scope of the present paper, and is left for a future publication), proves that, next to the spin glass 
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transition, q is a function independent of temperature, and thus nontrivial even at criticality. 
The prefactor q\, however, disappears at T c ensuring continuity of the order parameter through 
the spin glass transition. 



We are extremely grateful to Imre Kondor for his thorough review of the paper prior to publi- 
cation, and also for his useful suggestions. 
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